A NEW LOOK AT SHORT-TERM IMPLIED VOLATILITY IN ASSET 

PRICE MODELS WITH JUMPS 



ALEKSANDAR MIJATOVIC AND PETER TANKOV 

Abstract. We analyse the behaviour of the implied volatility smile for options close to expiry in 
the exponential Levy class of asset price models with jumps. We introduce a new renormalisation 
of the strike variable with the property that the implied volatility converges to a non-constant 
limiting shape, which is a function of both the diffusion component of the process and the jump 
activity (Blumenthal-Getoor) index of the jump component. Our limiting implied volatility 
formula relates the jump activity of the underlying asset price process to the short end of 
the implied volatility surface and sheds new light on the difference between finite and infinite 
variation jumps from the viewpoint of option prices: in the latter, the wings of the limiting 
smile are determined by the jump activity indices of the positive and negative jumps, whereas 
in the former, the wings have a constant model-independent slope. This result gives a theoretical 
justification for the preference of the infinite variation Levy models over the finite variation ones 
in the calibration based on short-maturity option prices. 

1. Introduction 

In financial markets, the price of a vanilla call or put option on a risky asset with strike e k and 
maturity t is often quoted in terms of the implied volatility a(t, k) (see 
definition and |10| for more information on implied volatility). Similarly, given a risk- neutral 
pricing model, one can define a function (t, k) \— > a(t, k) via the prices of the vanilla options 
under that model. The implied volatility is a central object in option markets and it is therefore 
not surprising that understanding the properties and computing the function (t,k) t-t a(t,k) 
for widely used pricing models has been of considerable interest in the mathematical finance 
literature. Typically, for a given modelling framework, the implied volatility a(t, k) is not 
available in closed form. Hence the study of the asymptotic behaviour in a variety of asymptotic 
regimes (e.g. fixed t and k — > ±00 [HI El [11]; t — > 00 with k constant [22] or proportional [13] 
to t; t — > and k constant [18j 2 1 . 7j etc.) has attracted a lot of attention in the recent years. 

In this paper we assume that the returns of the risky asset S = e are modelled by a Levy 
process X and study the relationship between the jump activity of X and the implied volatility 
at short maturities in the model S. Most existing approaches analyse either the at-the-money 
case, when the implied volatility is determined exclusively by the diffusion component and 



12) in Section for the 
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FIGURE 1. The liquid 10-A and 25-A strikes (left panel) and the corresponding implied 
volatilities (right panel) for market defined maturities t 6 {1 day, 1 week, 2 weeks, 1 month} 
in USDJPY suggest the following: as maturity t becomes small, the relevant strikes k t should 
approach the at-the-money strike and the implied volatilities a(t, kt) should remain bounded. 



converges to zero in the pure jump models (see [211 Prop. 5], [EJ [12]), or the fixed-strike out- 
of-the-money case, when the implied volatility for short maturities explodes in the presence 
of jumps (|17|. [7], [H] ). However, in the option markets, (a) although the implied volatility 
for liquid strikes grows with decreasing t, it remains within a range of reasonable values and 
appears not to explode, and (b) the liquid strikes become concentrated around the money 
as the maturity gets shorter. For instance, in the FX option markets, which are among the 
most liquid derivatives markets in the world, options with fixed values of the Black-Scholes 
delta are quoted for each maturity (see [2] for the details on the conventions in FX option 
markets and a natural parameterisation of the smile using the Black-Scholes delta). The market 
data in Figure [T] therefore suggests that, in order to understand the behaviour of the volatility 
surface at short maturities, one should look for a moving log-strike kt ^ 0, for t > 0, such that 
(i) the corresponding implied volatility has a non-trivial limit limq.n &(t, kt) and (ii) the strike 
kt converges to the at-the-money strike as maturity t tends to zero (i.e. lim^o kt = 0). 

This paper defines a new universal and model-free parameterisation of the log-strike given by 

k t = Byjt log(l/t) where 9 6 M\{0}. 

For fixed 9, the corresponding strike value tends to the at-the-money strike as t J, but is 
out-of-the- money for each short maturity t > 0. We prove that under suitable assumptions the 
limiting implied volatility c (f5) = lim^o &(t, h) takes the following form as a function of 9: 

(1) o- (9) = max I - ~ e ,a, — f =i for any 9 G M\{0}. 
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In this formula a denotes the volatility of the Gaussian component of the underlying Levy process 
X and a+ (resp. a_) denotes the jump activity (Blumenthal-Getoor) index of the positive (resp. 
negative) jumps of X. More precisely, if the jump measure of X is denoted by a + and a_ 
are given by 



Unlike in the case of fixed strike, where short maturity smile explodes in the presence of 
jumps, our parameterisation of the strike as a function of time yields a non-constant formula for 
the limiting implied volatility, which depends on the balance between the size of the Gaussian 
volatility parameter and the activity of small jumps. It allows us to make the following obser- 
vations about the relationship between the short-dated option prices and the characteristics of 
the underlying model: 

(i) the formula for (Tq{9) depends on the jump measure of the log-spot process X only if the 
jumps are of infinite variation; put differently, if the jumps of X are of finite variation, 
then the absolute value of the slope of the limiting smile for large \9\ is equal to one and 
in particular ao(8) does not depend on the structure of jumps; 

(ii) the limiting smile <jq(9) is V-shaped in the absence of the diffusion component (i.e. when 
(7 = 0) and is U-shaped otherwise; 

Remark (jij) provides a theoretical basis for distinguishing between the models with jumps of 
finite and infinite variation in terms of the observed prices of vanilla options with short maturity. 
It is well-known that, for any short maturity t, the market implied smile k \— > a(t,k) exhibits 
pronounced skewness and/or curvature, due, in particular, to the risk of large moves over short 
time horizons perceived by the investors. Hence, jumps are typically introduced into the risk- 
neutral pricing models with the aim to capture this risk and modulate the at-the-money skew 
of the implied volatility a(t, k) at small t (see e.g. [101 Eq (5.10)]). However, since this task can 
be accomplished by jumps of either finite or infinite variation, this requirement tells us little 
about the options implied jump activity of the underlying risk-neutral model. On the other 
hand, the formula for ctq{6) implies that, if we need to control the tails (in the parameter 9) of 
the implied volatility for short maturities, we must use jumps of infinite variation. This finding 
complements the analysis in [6] of the path- wise structure of the risk-neutral process implied by 
the option prices on the S&P 500 index. 

In recent years, there has been a lot of interest in the literature on the statistics of stochastic 
process in the question of the estimation of the Blumenthal-Getoor index of models with jumps 
based on high-frequency data. For example, it is shown in pQ that the jump activity (measured 
by the Blumenthal-Getoor index) estimated on high-frequency stock returns for two large US 
corporates is well beyond one, implying that the underlying model for stock returns should have 




and 
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jumps of infinite variation. Likewise, the formula in ([I]) suggests that jumps of infinite variation 
are needed in order to capture the correct tails (in 0) of the quoted short-dated option prices. 

The formula in ([!]) follows from Corollary |4j which gives the expansion of the implied volatility 
a(t,kt), where k t = 6y/t log(l/i), up to order o (1/ log(l/i)). This expansion is consequence of 
(A) Theorem [3j which itself gives an expansion of the implied volatility for a general log-strike 
kt that tends to zero as t \, 0, and (B) Theorem [I] and Proposition [2j which describe the 
asymptotic behaviour of the option prices under Levy processes with infinite and finite jump 
variations respectively. Theorem [3] relates the asymptotic behaviour of the vanilla option prices 
under a general semimartingale model to the asymptotic behaviour of the implied volatility as 
the log-strike kt tends to zero (it should be noted that the asymptotic regime (t, kt) in Theorem|3] 
is not covered by the analysis in [9], see Remark ( |rv| ) after Theorem [3] for more details). The 
asymptotic formula in Corollary [4] then follows by combining Theorem [3] with the asymptotic 
behaviour of the vanilla option prices established in Theorem [T] (for the case of jumps of infinite 
variation) and Proposition [2] (for jumps of finite variation). 

In a certain sense, Theorem [T] and Proposition [2] represent the main contributions of this 
paper. The asymptotic formulae for the call and put options, struck at e kt and e~ kt respectively, 
have the same structure in both results: the leading order term is a sum of two contributions, 
one coming from the diffusion component of the process and the other from the jump measure. 
Which of the two summands dominates in the limit depends on the level of the parameter 9. 
This structure of the asymptotic formulae is also reflected in the expression for o"o(#), as it is 
clear from ([l]) that <7q{6) = a if 9 is between — cryl — (a_ — 1) + and (J\J\ — (a + — 1) + , and 
&o(d) only depends on the jump measure otherwise. However, the proofs of Theorem [I] and 
Proposition [2] differ greatly: the finite variation case follows from the Ito-Tanaka formula, which 
can in this case be applied directly to the hockey-stick payoff function, while the case of jumps 
with infinite variation requires a detailed analysis of the asymptotic behaviour of the option 
prices. 

The remainder of the paper is organised as follows. Section [2] defines the setting and states 
Theorem [T] and Proposition [2] In Section [3j we state and prove the asymptotic formulae for 
the implied volatility and establish the limit in ([!]). Section [4] presents numerical results that 
demonstrate the convergence of option prices and implied volatilities given in the previous two 
sections, in the context of a CGMY model and a CGMY model with an additional diffusion 
component. Section [5] concludes the paper by proving Theorem [T] and Proposition [2} The 
appendix contains a short technical lemma, which is applied in Section [5] 

2. Option price asymptotics close to the money 



In this paper we study the behaviour of option prices close to maturity in an exponential Levy 
model S = e , where X is a Levy process with the characteristic triplet (a 2 , ^,7). Throughout 
the paper we assume the following: 
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• S is a true martingale (i.e. the interest rates and dividend yields are equal); 

• S is normalised to start at So = 1 (i.e. as usual the Levy process X starts at Xq = 0); 

• the tails of the Levy measure v admit exponential moments: 

(2) / e\ z \ {l+8) u{dz) < oo for some 5 > 0. 

J\z\>l 

In particular, assumption ^ guarantees the finiteness of vanilla option prices for any maturity 



t > 0. Section 2.1 describes the asymptotic behaviour of option prices for short maturities in 



the case the process X has jumps of infinite variation. Section 2.2 deals with the case where the 
pure-jump part of X has finite variation. 

2.1. Levy processes with jumps of infinite variation. Theorem [l] describes the asymp- 
totic behaviour of option prices in the case the tails of the Levy measure of X around zero 
have asymptotic power-like behaviour. This assumption does not exclude any exponential Levy 
models that appear in the literature but yields sufficient analytical tractability to characterise a 
non-trivial limit as maturity tends to zero for the option prices around the at-the-money. Before 
stating the theorem, we recall standard notation used throughout the paper: functions f(t) and 
g(t), where git) > for all small t > 0, satisfy 

fit) 

(3a) fit) ~ git) as 1 1 if lim ^ = 1, 

40 git) 

fit) 

(3b) fit) = oigit)) as 1 1 if lim ^ = 0, 

40 git) 

fit) 

(3c) f(t) = 0(g(t)) as 1 1 if is bounded for all small t > 0. 

Furthermore we denote x + := max{x,0} for any i£i 

Theorem 1. Let X be a Levy process as described at the beginning of the section and assume 
that the following holds 

(4) \\mx a+ v((x, oo)) = c_|_, \\m.x a ~i>({— oo, — x)) = c_ 

for a + ,a_ € (1,2) andc+,C- G [0, oo). Let kt be a deterministic function satisfying 

k t >0 Vt > 0, limfct = 
40 



and 



t l/a 

if a 2 = 0, lim— — = for some a G (max(a_, a + ), 2), 
40 kt 

ifa 2 >0, lim^=0. 
40 k t 
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Then, if c + > 0, we have 

(5) E[(e x * - e kt )+] ~ E[(e CT ^- V - e *t)+] + ^ °± as t ; 

a + — 1 

and, if C- > 0, it holds 

(6) E [(e- kt -e Xt ) + ]^E[{e~ kt -e' 7Wt -^) + } + t ^ ^ as 1 10. 

a_ — 1 

Remarks. (i) Theorem [I] implies that the price of a call (resp. put) option struck at e kt (resp. 
e~ kt ) tends to zero at a rate strictly slower than t if the paths of the pure jump part of X 



have infinite variation. In particular, combining the notation in (3a) and (3b), we get that 
the following equalities hold as t \. 0: 

E[(e Xt - e kt ) + ] = E[(e CTiy '-V - e kt ) + ] + -± + o(tk]- a+ 

a+ — 1 V 

l-a- 



E [( e -*t - e x *)+] = E[(e- fct - e aWt ~^) + ) + ^ ~ C ~ +o(tk\- a - 



(ii) The proof of Theorem [T] is given in Section 5.1 



2.1.1. Blumenthal-Getoor index and the short-dated option prices. Recall that for any Levy 
process Y with a non-trivial Levy measure iy, the Blumenthal-Getoor index, introduced in [3], 
is defined as 



BG(Y) := inf J r > : / 

1 



(7) BG(F) := inf <( r > : / < oo } . 

(-i,i)\{o> J 

The Blumenthal-Getoor index is a measure of the jump activity of the Levy process Y, since the 
following holds: r > BG(Y) if and only if ^2 s<t \AY s \ r < oo almost surely, where AY, := Y s — Y s - 
denotes the size of the jump of Y at time s. Furthermore, it is clear from (J7| that BG(y) lies 
in the interval [0, 2]. 

In recent years there has been renewed interest in the Blumenthal-Getoor index from the 
point of view of estimation of the jump structure of stochastic processes based on high-frequency 
financial data. For example, it was estimated in pQ that the value of BG(Y) is around 1.7 (i.e. 
the stock price process has jumps of infinite variation) based on high-frequency transactions 
(taken at 5 and 15 time intervals) for Intel and Microsoft stocks throughout 2006. 

Let X + and X~ be the pure-jump parts of the Levy process X from Theorem [l] In other 
words X + (resp. X~) is a Levy process with the characteristic triplet (0, v + , 0) (resp. (0, v~ , 0)), 
where u + {dx) := l{ x>0 yi/(dx) (resp. v~{dx) := l{ x< Q}v(dx)). Then assumption Q implies 

BG(X+) = a+ and BG(X~) = a_, 

and relations ^ and ^ of Theorem [l] describe how the Blumenthal-Getoor indices of the 
positive and negative jumps of X influence the asymptotic behaviour of option prices at short 
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maturities. The result clearly depends on the asymptotic behaviour of the log-strike kt- In 
Section [3] we will prescribe a specific parametric form of kt (see (13)) and give explicit formulae 



for the asymptotic expansion and the limit of the implied volatility as maturity tends to zero in 
terms of the Blumenthal-Getoor indices of X + and X~ (see Corollary [4] for details). 

2.2. Levy processes with jumps of finite variation. In this section we study the option 
price asymptotics at short maturities in the case the process X has a (possibly trivial) Brownian 
component and a pure jump part of finite variation. 

Proposition 2. Let X be a Levy process as described at the beginning of Section Assume 
further that the jump part of X has finite variation, i.e. 



\x\v(dx) < 00. 
Let kt be a deterministic function satisfying 



and 



k t >0 Vt > 0, lim k t 
40 



ifa 2 = 0, lim— = 0, 



ifa 2 >0, lim^ = 0. 



40 kt 
40 kt 



Then, as t \, 0, it holds: 



(8) E[(e Xt - e kt ) + ] = E[(e aWt ~ V - e fe *)+] + t I {e x - l)v(dx) + o(t) 

J (0,0a) 

and 

(9) E[(e~ kt - e Xt ) + ] = E[(e~ kt - f f w *- s r)+] +t f (1 - e x )u{dx) + o{t). 

Remarks. (i) Proposition [2] implies that, in the absence of a Brownian component, the call 
and put prices of options struck at e kt and e~ kt , respectively, tend to zero at the rate equal 
to t if X has paths of finite variation (cf. Remark ^ after Theorem [I]). 
(ii) The Blumenthal-Getoor indices of the positive and negative jump processes X + and X~ 



of X, defined in Section 2.1.1 are both smaller or equal to one by the assumption in 
Proposition [2] Furthermore, unlike in the case of jumps of infinite variation, Proposition [2] 
implies that the asymptotic behaviour of short-dated option prices (as maturity t tends to 
zero) does not depend up to order o(t) on the indices BG(X + ) and BG(X~). Hence, the 
same will hold for the short-dated implied volatility (cf. Corollary El) . 



(iii) It should be stressed that the proof of Proposition [2j given in Section 5.2, is fundamentally 



different from that of Theorem [TJ as it relies on the path- wise version of the Ito-Tanaka 
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formula for the processes of finite variation, which cannot be applied in the context of 
Theorem [TJ 

3. Asymptotic behaviour of implied volatility 

The value C BS (t, ft, a) of the European call option with strike e k (for any k 6 M) and expiry 
t under a Black-Scholes model (with log-spot Xt = crWt — ta 2 /2 of constant volatility a > 0) is 
given by the Black-Scholes formula 

(10) C BS (t,k,a) = N(d + )-e k N(d_), where d± = 

(Jyjt 2 

and iV(-) is the standard normal cumulative distribution function. The price of a put option 
with the same strike and maturity is given by P (t, k, a) = e k N(—d+) — N(-d-). Let S be a 
positive martingale, with Sq = 1, that models a risky security and denote by 



;il) C(t,k):=E 



S t -e k 



and P(t, k) := E 



e k -S i 



the prices of call and put options on S struck at e k with maturity t, respectively. The implied 
volatility in the model S for any log-strike k G IR and maturity t > is the unique positive 
number a(t, k) that satisfies the following equation in a: 

(12) C BS (t,k,a) = C(t,k). 

Implied volatility is well-defined since the function a i— > C BS (t,k,a) is strictly increasing on the 



positive half-line and the right-hand side of (12) lies in the image of the Black-Scholes formula 



by a simple no- arbitrage argument. Put-call parity, which holds since S is a true martingale, 
implies the identity P BS (t, k,a(t, ft)) = P(t, k). 

In order to study the limiting behaviour of the implied volatility close to the at-the-money 
strike 1 = e° for short maturities, we define the following parameterisation of the log-strike kt- 



1 



1/2 



(13) kt := e\t\og-\ , where 9 € M\{0}. 

We can now define the implied volatility at '■ K\{0} — > (0, oo) as a function of 9 in the asymptotic 



maturity-strike regime (t, ftt), given by (13), for a short maturity t: 



(14) <T t (0) := a(t,k 



The implied volatility crt(9) is of interest in the context of processes with jumps, because its 
limit <tq{9), as 1 1 0, exists and is finite for each 9, depends on both the jump and the diffusion 
components of the process and can be computed explicitly in terms of the parameters. In order 
to find the asymptotic behaviour of o~t(9), we first state Theorem[3j which relates the asymptotics 
of <Jt{9) to the asymptotic behaviour of the out-of-the- money option price 

(15) ItiO) := C(t, h)l {e>0} + P(t, h)l {e<0} 
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under the model S as maturity t tends to zero. 

Theorem 3. Let S be a martingale model for a risky security with Sq = 1 and kt a log- 



strike given in (13) for a fixed 9 £ M\{0}. Let Ct and Pt be deterministic functions such that 
C(t, kt) ~ Ct and P(t, kt) ~ Pt as t \. 0, where C(t, kt) and P(t, kt) are given in ( 11 ), and define 
lt{0) := Cil|g >0 } +i 3 fl{0 < o}- Assume further that the out-of-the-money option price lt(0), given 
in (15), satisfies: 

(lb) - < nmml < limsup < oo. 

2 40 log* 40 logt 



Then the implied volatility <Jt(6), defined in (14), can be expressed by 



(17) 
and 

(18) 



\e\ |0|iog (2£ 'V' V27r 



y/2L t (0) - 1 (2L t {d) 



log} 



+ o 



log 



21 1 



+ 



1^1 lo s w — 



2L t (9) - 1 



2L t {9) 



+ o 



lOg i I log J 



where Lt{9) := Jt(lt(0)) and Lt(8) := Jt(lt{0)) are defined by the formula 



Jt{x) :-- 



log x log log \ 



for any x, t > 0. 



as 1 1 0, 



as 1 10, 



log t log \ 

Before proceeding with the application and proof of Theorem [3j we make the following remarks 
in order to place it in context. 

Remarks. (i) In the Black-Scholes model with volatility a > 0, the following well-known ex- 
pansion of the call option price in the (t, kt) maturity-strike regime (13) holds (e.g. a 
straightforward calculation using [9, Eq. (3.10)] yields the expansion): 



(19) C BS (t,k t ,a) 



: £2 + 2a 



1 



1 



'2ir 



e 2 lop 



e 4 ios 



t + o 



i 



log 3 I 



as 1 10. 



In particular we have log C BS (t, k t ,a) = (| + Jf?) logt + o(log(l/t)) as t J, and hence the 
assumption in (16) is satisfied in the Black-Scholes model, 
(ii) Note that the log-strike kt in (13) satisfies the assumptions of Theorem [TJ For any Levy 
process X as in Theorem [TJ formula ^ and Remark ^ above imply 



(20) 



log C(t, kt) = min 



a+ 1 

2 ' 2 2ct 2 j 



1 



logt + o(log(l/i)) as tiO. 



Since the minimum of the constants in front of logt is clearly larger than 1/2, assump- 
tion (16) of Theorem [3] is satisfied. As we shall soon see, it is the balance (as a function of 
9) between the two constants in (20) that determines the value of the limiting smile ctq(0). 
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(iii) Let a Levy process X be as in Proposition [2] (i.e. with jumps of finite variation). For- 
mulae ((8| and (19) imply that the call option price C(t,k t ) under the model S has the 



following asymptotic behaviour 



(21) logC(t,k t 



min si, — I « 

2 2<r 2 



logt + o(log(l/i)) as iiO. 



In particular note that assumption (16) is satisfied and that, in the case of jumps with 
finite variation, the constant in front of logt does not depend on the Levy measure but 
solely on the diffusion component of the model, 
(iv) In [9j the authors present a general result, which translates the asymptotic behaviour of the 
option prices, in a generic maturity-strike regime, to the asymptotics of the corresponding 
implied volatilities. Unfortunately the results in [9] do not apply in the regime (t, hi), for kt 
in (13), since the standing assumption of [9], max{0, log(l//^)} = o(log(l/C(t, kt))) (see [H 
Eq. (4.3)]), is not satisfied in our setting by (20) and plj ). We therefore have to establish 
Theorem [3j which is applicable in our context as remarked in (JTTJ) and (iii) above. 

Before proving Theorem [3j we apply it, together with Theorem [l] and Proposition [2j to derive 
the main asymptotic formula of the paper. 

Corollary 4. Let X be a Levy process with the jump measure v and the Gaussian component 
a 2 > 0. Pick 9 G K\ {0} ; let kt be the log-strike from (13) and let o~t(0) be the implied volatility 
defined in ( 14 ) . Then the following statements hold. 

(a) Let X be a Levy process satisfying the assumptions of Theorem^ Then the implied volatility 
o~t{0) takes the form 

(22) 

±d [l + I ± (t,9)] + o 




if ±9 > o-\J2 — ct± and c± > 0, 
if < ±9 < a^2 - a± and c± > 0, 



3 - a± log log g 
2(2 -a±) log} 




a±) 2 c±v27r \ 1 



log 



1 - 
t 



as t ],0, 



for t > 0, 



and the sign ± denotes either + or — throughout the formulae in (22) and (23). In par- 
ticular, the limiting smile oq{9) := lim^o &t{6) exists for any 6 £ 1\ {0} and takes the 
form 



M0) 



max 



±9 



if c± > 0. 



V2 - «± ' 

(b) Let a Levy process X be as in Proposition^ and let 7+,7_ > be equal to the following 
integrals 



7+ 



(0,00) 



l)v(dx) 



7- 



(-00,0) 



(1 -e x )v{dx). 
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Then the implied volatility o~t(9) for short maturity t is given by 



(24) a t {9) 
where 

(25) 



±9 [l + F±(t,0)]+o K ^r ,. 



a + o 



T I > 

t 



if ±9 > a and y± > 0, 

if < ±9 < a and y± > 0, 



as 1 10, 



F±(t,9) := — j-t. +log 

lo gl 



7-tV27r \ 1 



/or t > 0, 



and ± denotes either + or — throughout the formulae in (24) and (25). The limit of the 
implied volatility smile as maturity tends to zero, ctq{9) := lim^o °~t(9), exists for 9 6 M.\ {0} 
and is egua/ to 

ao(9) = max{±9,a} if y± > 0. 

Remarks. (i) Recall display Q in Theorem [l] and note that the assumptions c+ > and 
c_ > of Corollary [4] Q mean that, as x J, 0, the tails around zero of the Levy measure v 
of X behave as v((x, 00)) ~ c+x~ a+ and 00, — x)) ~ c_:r~ a ~. Note further that, once 
we have identified the precise rate of the tail behaviour of v at zero, the constants c+ and 
c_ do not feature in the limiting formula o~o(9). 
(ii) The assumption y± > in Corollary [4] Q ensures that the process X has positive jumps 
when 9 > and negative jumps when < as we are only interested in the asymptotic 
behaviour of the implied volatility in the presence of jumps. 

In Section 3.1 we derive Corollary [4] from Theorem[3]and in Section 3.2 we establish Theorem[3j 



3.1. Proof of Corollary [3[ (a) Assume first that o-y/2 - a + > 9 > 0. Define C t 
C BS (i, kt, o~) and note that ( |19[ ), the definition of kt in (13) and ^ of Theorem [l] imply 



(26) C(t, k t ) ~ C t) 



and hence 



log C(t,fct) log C t 



logt 



logt 



+ 



1 



log I 



as t 4 0, 



where C(t,kt) denotes the call option price with maturity t and strike e kt under the exponen- 
tial Levy model e x . Assumption (16) of Theorem [3] is therefore satisfied by Remark (|I| after 
Theorem 3l The formula for Lt{9) = log Ct/ logt — (log log V) log | takes the form 



(27) 



Lt{9) 



2 + 2a 2 



log 



(7 



9 2 V2vJ log 7 



+ 



log 



1 ' 



as 1 1 0, 



The formula in (18) of Theorem [3j together with (27) and the Taylor expansions in log(l/i) as 

no 



2L t {9) - 1 



1 



1 °\ ( ^ 
1 + F l0g Uvir7 logf 



+ 



log 



1 ' 
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9 log 



(2£t(g)-l) 



3 



a 3 log 



9W2n 



2L t (9) - 1 



log I 



9 2 



+ o 



log 1 \ log 7 



yield the formula in ( 22 ) . 



In the case — a + < 9, the relation (|26|) is satisfied by Ct ■ 



. This follows directly 

from the definition of kt in (13) and Theorem [T] (see formula ([5])). An analogous argument as the 
one above shows that in this case the assumptions of Theorem [3] are also satisfied. By definition 
of Lt(6) in Theorem [3J we find 

3-a+loglogj : 



2Lt(0) - 1 = (2 - a+) 



1 



log 



jl-«4 



2 — a+ log j 2 — a + 
By Taylor's formula the following asymptotic relations hold as 1 4- 0: 



a + — 1 / log 



1 



'2L t {9) 
and 



y/2-a+ 



1 + 



2(2 



log log! 

log^ 



+ 



1 



a j 



log 



9 i- a+c+ 
a+ — 1 



1 

log! 



+ 




01og 



(2L t (0)-l 



3 



01og 



(2-a + )gV2^ 



2Lt(0) - 1 



log 



(2 



a+) 2 



+ 



lOg I llOgi 



Substituting these expressions into (18) establishes the formula in ( |22[ ), 

Assume now that —o^J2 — a_ < < 0. Define P t := P BS (t,k t ,a), where P BS (t,k t ,a) is the 
put option price in the Black-Scholes model, and recall the well-known put-call symmetry 

> BS (t,k t ,a) 



(28) 



P L 



^C BS (t,-k t ,a), 



which holds since the laws of minus the log-spot under the share measure (i.e. the pricing 
measure where the risky asset is a numeraire) and the log-spot under the risk-neutral measure 
(i.e. the measure where the riskless asset is the numeraire) coincide. Analogous to the case 
above, (19) with the put-call symmetry, the definition of kt in (13) and Q of Theorem [T] imply 



(29) P{t, k t ) ~ P t , and hence 



log P(t,k t ) _ log P t 



log* 



log* 



+ 



log^ 



as 1 1 0, 



where P(t,kt) is the put option price under the exponential Levy model e . Therefore the 
assumptions of Theorem [3] are satisfied and Lt(0) takes the form (27). Note that the right-hand 
side of (27) depends solely on the even powers of 9 and hence the fact 9 < does not influence 



the asymptotic behaviour of Lt{9). The proof of formula (22) now follows in the same way as 
in the call case above. 

-a^2-a- > 9 we define P t — <( ~ fct)1 "" +c+ 



In the case 



Q+ — 1 



Under this assumption, the rela- 



tion (29) is satisfied by ^ of Theorem [I] and the rest of the proof follows along the same lines 
as in the case Oy/2 — a+ < 9. This proves formula (22). 
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(b) The proof of part (b) of the corollary is based on Proposition [2] and Theorem [3j The steps 
are analogous to the ones in the proof of part (a): 

if a > > 0, define C t := C BS (i, k t ,a); if a < 9, define C t := 

if - a < 9 < 0, define P t := P BS (t, k t , a); if -a>9, define P t := fry_. 

The details of the calculations are left to the reader. 

3.2. Proof of Theorem [3| We first assume that 9 > 0. Equality (19) implies the following 
log C BS (t,k t ,s) loglogi_l 9 2 1 s 3 Q s 2 1 ( 1 ^ 



logt log I "2 2s 2 log I 6 9 2 V2^ 2 log 2 \ \log 3 \ 

as t J. for any s > 0. Define 

F(s,z) := -zlogC BS (e- 1/2 ,A: e _iA,s) + ^log^ 



and note that F(s, z) corresponds to the left-hand side of the above formula with the change 
of variable z = 
equality holds 



of variable z = j^p - - The expansion shows that F(s,z) is regular as z — > and the following 



The expansion for the inverse mapping can be deduced from this expression as follows. To keep 
the formulae simple, we give the expansion up to 0(z 2 ): 

1 9 2 s 3 
F(s, z) = a(s) + zb(s) + 0(z 2 ), where a(s) = - + b(s) = log 



2 2s 2 ' w 9 2 V2^' 
Denote by F~ l (y,z) the unique positive solution of the equation F(s,z) = y, where y equals 
J e -i/z(x) (see the statement of Theorem [3] for the definition of Jt{x)) and x is any arbitrage-free 
call option price with maturity e _1//z and strike k e -i/ z . The uniqueness of the quantity F~ 1 {y, z) 
is equivalent to the fact that the implied volatility is a well defined quantity. 
An approximate expression for y is given by 

y = a{F-\y, z)) + z^F'^y, z)) + 0(z 2 ) 

and hence we find 

a-\y) = a-\a(F-\y, z)) + zb(F~ 1 (y, z)) + 0{z 2 )). 

Using the regularity of the coefficient a in the neighbourhood of the point F _1 (y, z) > 0, we can 
expand the inverse a -1 around the point a(F^ 1 (y, z)) as follows: 

a-\y) = F- l (y, z) + [pT x )' [a{F~\y, z)))b(F~\y, z))z + 0(z 2 ). 

In view of this expression, and using once again the regularity of the coefficients a and b, we can 
replace F^ 1 (y, z) with a~ l (y) in the second term, obtaining 

a-\y) = F-\y, z) + {a' 1 )' \y)b{ar\y))z + 0(z 2 ). 
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Hence, the following asymptotic equalities hold true: 



F-\y,z) = a-\y)- l *f^z + 0^) 
a' [a %)) 

a fling &t^)l^E 



(2y-i) 



Substituting the expression for F, we find an expansion for the implied volatility (Jt(9) given 
in ([17]). Now, C{t, k t ) ~ C t implies that log ^' fct) - = o(log- 1 r 1 ). Since all the coefficients 



log t log t 

in expansion ( |17[ ) are regular, the additional term arising from this difference may be ignored in 
an expansion up to order o(log _1 i _1 ) and (18) follows. 



The formulae in the theorem in the case 9 < will be established by applying the result 
for the positive log-strike under the share measure. More precisely, let P denote the original 
risk-neutral measure under which the process S is a positive martingale started at one. For each 
time t, we define the share measure P on the cr-algebra Tt of events that can occur up to time t 
via its Radon-Nikodym derivative 3p|j" t ■= St and note that the following relationship holds for 
any log-strike k € M: 



(30) 



P(t,k) =E 



St) 



(Si 1 



)+ = e k C (t,-k), 



where C(t, —k) := E \(S t 1 — e k ) + ] denotes the expectation under the share measure P of a call 
payoff with strike e~ k , where the evolution of the risky asset is given by S' 1 . Note that S' 1 is a 
positive martingale started at one under P and hence C(t, —k) represents and arbitrage-free call 
option price. Furthermore, the put-call symmetry formula in the Black-Scholes model (see (|28[)) 



and the equality in (30) mean that the implied volatility a(t, k) defined by the put price P(t, k) 
coincides with the implied volatility a(t, —k) defined by the call price C(t, —k) (see beginning 
of Section [3] for the definition of a(t, k)). 

Note that, since 9 < 0, we now have —kt > and at(9) = a t (—9), where denotes 
cr(i, —kt). In order to apply the formula in (17) to C(t, —kt), we have to ensure that assump- 



tion (16) is satisfied. Since (16) holds for P(t,k t ) and kt = o(logi), the equality in (30) im- 
plies (16) for C(t, —k t ). Therefore formula (fl7|) gives an asymptotic expansion of <Jt(9) = at(—9) 



in terms of Lt (- 



Jt(C(t, —kt)). Since equality (30) implies 



L t {9) = L t (-9) - 6y/t/ log(l/t) = L t (-9) + O 



log 



2 1 



as 1 1 



and the two leading order terms in (17) are regular in Lt(—6), the asymptotic expansion in (17) 
also holds when Lt(—9) is replaced by Lt(9). The formula in (18) now follows by the same 
argument as in the case of the positive log-strike. This concludes the proof of the theorem. 
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4. Numerical results 



In this section, we present some numerical illustrations for the convergence results discussed 
in Section [3j We focus on the generalised tempered stable Levy process X with Levy density 

-\+x „ „— A— \x\ 



(31) 



v(x) 



c_|_e 



|i+a+ 1 {^>o} 



+ 



c_e 



|l+a- 



-1 



{x<0}- 



This class of processes includes the widely used CGMY models (see e.g. [1]). For this process, 
the price of a European call option with pay-off (Soe Xt — K) + at time t can be computed as 

iu—R 

(32) 



K 
2^ 



<f>t(-u - iR) 



-du. 



(R - iu)(R - 1 - m) 

where <f>t is the characteristic function of X t and R > 1 (see e.g. [5] or [21]). We compute the 



integral in ( 32 ) with an adaptive integration algorithm. 



4.1. Testing the algorithm. To ensure that the prices returned by our algorithm are 
correct, we first compare them to the values computed in [23J with their approximate 
"fixed point" algorithm (PDE discretisation). The following table shows that the values we 
obtain are very similar with the small discrepancy probably due to the discretisation error of [23J . 



s 


K 


T 


r 


c = c + = c_ 


A+ 


A_ 


a = q+ = a_ 


Value ([23]) 


Our value 


90 


98 


0.25 


0.06 


16.97 


29.97 


7.08 


0.6442 


16.212578 


16.211904 


90 


98 


0.25 


0.06 


0.42 


191.2 


4.37 


1.0102 


2.2307031 


2.2306558 


10 


10 


0.25 


0.1 


1 


9.2 


8.8 


1.8 


4.3714972 


4.3898433 



4.2. Convergence of the at-the-money (ATM) options. In this section we fix the param- 
eters of the tempered stable process at 



(33) 



c+ = c_ = c = 1, A + = A_ = 3, a+ = a_ = a = 1.5 



and So = 1. First we analyse the rate of convergence to zero of the ATM options. It follows 
from the results in [15] that the ATM option price satisfies 

E[(e Xt - 1) + ] ~ t 1/a E[(Z*) + ], 

where Z* is a stable random variable with the Levy density , . Furthermore it is known that 

\x\ 



E[(Z* 



(2c) 



l/a 



, „, . ira\ 1 / a 
r(l - l/a) -r(a)cos — I 



=: C. 



it ■ ■ \ • • 2 J 

Figure 2 plots the dependence of the normalised option price r 1 / Q E[(e Xt - 1)+] and the nor- 
malised "Bachelier" price t- 1 / Q E[(X<) + ] on logt, i.e. on time to maturity expressed on the 
log-scale. The horizontal line in Figure [2] corresponds to the value of the constant C. The 
desired convergence is clearly visible. 
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Figure 2. Convergence of the re- normalised price for ATM options: the param- 
eter values of the tempered stable process, which models the log-stock are given 
in (ph. 



4.3. Convergence of option prices with variable strike. In this section we investigate 
numerically the convergence of the out-of-the- money (OTM) option prices given in Theorem [T] 
The parameter values for the underlying process are given in ( 33 ) . Note that in the case of the 
tempered stable Levy process with Levy density (31), the limits in Q of Theorem [T] take the 
form 

c 
a 

Figure [3] shows the dependence of the normalised option and "Bachelier" prices, respectively 
given by 



limx a u((x, oo)) = limcc c V((— oo, — x)) 

x\.0 x\,0 



E[(e Xt -e kt Y 



on time to maturity in log-scale, where 

kt = t^' 



and 



with 



mx t - k t ) 



l-a 



a 



tk 



1.9. 



predicted by Theorem jlj 



The horizontal dotted line shows the limiting value , c ^ : 

° a(a—l) 

Similarly, Figure [4] plots the dependence on time to maturity (on the log-scale) of the nor- 
malised option price 



E[(e 



tk 



l-a 



for kt 



tlog 



and 



E {0.1,0.2,0.3}. 



As in Figure 3, the limiting horizontal dotted line is given be 



a(a— 1) 
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1.0- 



0.5- 



o.o 




: Normalized option price 
Normalized Bachelier price 



— Limit 



■10 



log(Time) 



Figure 3. Convergence of the re-normalised price for OTM options: kt 



with a' = 1.9 and the other parameters of the process X are given in (33). 




4.4. Convergence of the implied volatilities to the limiting smile. In this section we 
illustrate the convergence of the implied volatility (expressed as function of the re-normalised 
strike 6) to the limit <jq(6) given in Corollary [4j In order to test the formula both with and 
without the diffusion component we fix two models: the first is a pure jump tempered stable 
Levy process with the following parameter values 



c + = c_ = c = 0.01, A + = A_ = 3, a + = a_ = a = 1.5, 
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which correspond to the unit annualised volatility of about 14%. The second model is the same 
tempered stable process with an added diffusion component of volatility a = 0.2. 

Recall that the limiting formula for positive 6 is gq{0) = max{<7, Figure JHJ plots the 

right wing of the implied volatility smile (as function of 9) for different times to maturity when 
a diffusion component is present (left graph) and diffusion component is absent (right graph), 
together with the limiting shape ctq(6). The convergence to the limit is visible in both graphs but 
slow, because the error terms in Corollary [4] are logarithmic in time. Nevertheless, the following 
observations can be made already at "not such small" times: 

• The smile is remarkably stable in time, when it is expressed as function of the re- 
normalised variable 9. In particular, the slope of the wings predicted by Corollary [4] is 
achieved rather quickly. 

• The distinction between the U-shaped smile in the presence of a diffusion component 
and the V-shaped smile in the pure jump case, is clearly visible. 



4.5. Approximation of the implied volatility for small times to maturity. In this sec- 
tion we illustrate the approximation of the implied volatility at small times by the asymptotic 



formula (22). We take the same parameters of the tempered stable process as in Section 4.4 and 
consider the case a = (when the diffusion component is present, in the region where the pure 
jump component dominates, the asymptotic formula is the same, and in the diffusion-dominated 
region, there are no additional terms added to the constant limit). Figure [6] illustrates the quality 
of the approximation for t = 1 day and t = 0.1 days. 



5. Proofs 
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0.30- 

0.25- 

0.20- 

0.15 

0.10 

0.05 



0.00 



Implied volatility 

3-term assymptotic approximation 



0.40- 
0.35- 
0.30- 



Implied volatility 

3-term assymptotic approximation 




0.00 0.05 



0.10 0.15 0.20 

Theta 



0.25 0.30 



0.00H 

0.00 



0.05 0.10 0.15 0.20 0.25 0.30 

Theta 



Figure 6. Approximation of the implied volatilities by the asymptotic formula 



(22). Left: t = 1 day. Right: t = 0.1 days. 



5.1. Proof of Theorem [Tl By Lemma[5j to prove Theorem[TJ it is sufficient to show that 

a + — 1 



(34) E[(X t - kt)+] = E[(aW t - k t )+] + ^ °± + o{tk]- a+ + E[(aW t - k t ) + ]) 



as 1 4- for the call case and 

(35) E[(-h - X t )+) = E[(-k t - aW t ) + ] + ^ + o{tk]- a - + E[{-k t - aW t ) + ]) 

(\ — 1 



as t 4 for the put case. Note that (35 ) follows from (34) by a substitution X i— > —X. Therefore, 



from now on we concentrate on the proof of (34), assuming with no loss of generality that c+ > 0. 



Step 1. In this first step, we assume that v{{— oo,0)) = and would like to prove 
(36) E[(X t - k t )+] = E[(aW t - k t )+] + ^ °± + o{tk]- a+ ). 



Fix t > 0, and e > with e < let X t be a Levy process with no diffusion part, Levy measure 
i^(dx)l{o<x<efci} an d third component of the characteristic triplet 



7t = 7 - ; zv{dz) 

'(ekt,l] 

Let (C*)«>i De a sequence of i.i.d. random variables with probability distribution 

v{ dz ) l {z>ek t } 

u({z : z > eh})' 

and iV* a standard Poisson process with intensity := v{{z : z > ekt}). Furthermore we 
assume that X*, iV* and (C*)i>i are independent. Then the following equality in law holds 



(37) 



X t ±vWt + X* + Y,& 



8=1 



20 



ALEKSANDAR MIJATOVIC AND PETER TANKOV 



and it follows that 

(38) E [(X t - k t ) + ] = e" Xtt E [(aW t + X\ - k t ) 

(39) + \ t te~ Xtt E [(aW t + X\ + £ - k t ) + ] 



(40) 



k>2 



E 



aWt + Xt + ^l-h 



i=l 



As a preliminary computation, we deduce from the assumptions of the theorem that the following 

. definit 

z 2 v(dz) 



asymptotic behaviour holds as t J, (recall definition (3a)): 

2 , ■ , 2c 



(41) A t ~ c + (ek t )- a + , S t := 



(O.fcte] 



2 - a_ 



(42) 



(43) 



7t 



-, E[(X*) 2 ] = t 2 ( 7t ) 2 + tSt~^, E[g]~-^+-eA* 

- 1 - , — 



a+ — 1 



E[(X<) 4 ] = t / 
iff 



z V<fc) + 4t 2 7i 



4c+ 



/ z 3 z/(cte) + 3i 2 £ 2 + 6t 3 7 t 2 Si + i 4 7 4 

■/(0,fc t el 



4 - a+ 



To estimate the term in (38), we apply the argument inspired by Lemma 2 in [19]. In the 
current notation this implies 



(44) 



{ rh/t ) 
F[Xl > k t ] < exp l-t J r(z)dz \ , 



where r : [jt, oo) — >• M is the inverse function of s : [0, oo) — > M defined by 



s(x) = 74 + 
By Taylor's theorem, this function satisfies 



(0,M 



z(e za; - l)i/(dz). 



2ktex _ j 



fete 



This implies that 



and therefore, substituting this into (44) 
PLX| > fc t ] < exp j- 

(45) < exp \ - 



2(k t ef 

h - it t ( he 
2k t e 



log(l + s)ds 



log ^(W*-7t) 



From the assumptions of the theorem and (41)— (43), there exists ii > such that t < to implies 

i 

,2, 



W>* t ]S-p{-i.og(^!)H 



2etS/ 



< C(tA; t - Q+ )ft < C(tk; a+ f 
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for some constant C < oo. By similar arguments it can be shown that 

h 



(46) 



X\> 



Coming back to the estimation of (38), we first deal with the case a = 0. In this case, the 



Cauchy-Schwartz inequality allows to conclude that 



E 



(xt - h) 



< E 



{x t t ) 2 Ynx t t >k t ]^ = o{k t {tk- t 



because the first factor remains bounded by d4lh-(43). 



Let us now focus on the case a > 0. Let /(x) := f^°( z ~ x ) e 2 dz. The expectation in 

= ay/ffl. 



( 38 ) can be expressed as 



E[(aW t + X\ - k t 



f 



k t -Xi 



(7 



Vi 



By Taylor's formula, we then get 



E[(aW t + X\ - k t )+\ = ayftf (^) - f (^) E[Xf\ 



a 



E 



{x\Y (i - e)f 



h - exj 

aVi 



do 



E[(aW t - h) + ] + 7ttF[aW t > h] + 



1 



axThrt 



=E 



I? 



1 { K t-v^i 
tri/t 



de 







We now need to show that the second and the third terms do not contribute to the limit. Since 
by assumption ^ ->■ 0, we have that P[crWt > kt] — > as t ->■ 0, and therefore, by ([41])- ([43]), 

IttnaWt > k t } = o(tk^ a+ ). 



The last term can be split into two terms, which are easy to estimate using (|41|)-(|43|): 
1 



Vt 



E 



(X t ) l {xl< k ±} I C 1 



1 k t -ext 

2 \ o^Tt 



de 



i 



< -jM[{XlY\e 



if fc t 



8 \crVt 



0(tk 



Vi 



o(tk]- a+ ) 



because by assumption of the theorem, % — >• oo. On the other hand, 



Vi 



E 



1 / Ki-»^t 



de 



< ^E 
Vi 



< ^fUlYMxi > |]l = 0(*f ^)0((tC*+) 2 ) = (^-«+) 



by (43) and (46). We have therefore shown that 



E[{aW t + X\ - k t ) + ] = E[(aW t - k t ) + ] + o(tk 



1 — OLA 
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From (41), the assumption on kt in Theorem [I] and the Lipschitz property of the function 
i4x + , it follows that 



e- Xit E[(aW t + X\ - k t ) + ] = E[(aW t - k t ) + ] + o{tk' t 



1— OLA 



as well. 



For the term in (39), the Lipschitz property of the function i 4 i+, (41)-(43) and the 
assumption of the theorem (i.e. the first assumption on kt in Theorem [T] in the case a = and 
the second one otherwise) imply the following estimate: 



E 



(aw t + x\ + ei - h) 



E 



{ei-h) + <\tt{E[\xt\]+<Tn\wt\]} 



< Xtt{E[\X^ 2 + aVt} = O(Attisf) + crA t t§ = o(tk 



1— a+ \ 



as t — > 0. 



On the other hand, integration by parts implies 



XttE 



t 



(z - k t ) v(dz) 



t 



U(z)dz 



tk) a+ c_i_ 
- as t -> 0, 



fa 



where U(z) := v((z, oo)), which yields the second term in ([5]). 



a_i_ — 1 

To treat the summand in (40), observe that by (|41|~(|43|), for k > 2, 

(aW t + X\ +Y£~ hj <°Vi + E[|X|| 2 ] 1/2 + fcEKi] 
cjv 7 * + 0{t l l 2 k l - a+/2 ) + feO(fct) = jfcO(jfet). 



E 



and hence o{tk\ a+ ) 



Therefore, the summand in (40) is of order 0(kt\ 2 t 2 ) = 0(kt(tk' t ~ a+ ) 2 ) 

Step 2. We now treat the case when z/((— oo, 0)) ^ 0. Let X~ be a spectrally negative Levy 
process with zero mean and zero diffusion part and Y be a spectrally positive Levy process such 
that X - + y = X. Let /3 € (max(a + , a_), a) (where we take a = 2 is c > 0) and xt = t& . As 
before, we fix e > and let X* be a Levy process with no diffusion part, zero mean and Levy 
measure v(dx)l{_ ext < x< Q}, let jt = J(_oo -ext) zv {dz), let i be a sequence of i.i.d. random 
variables with probability distribution 

v{dz)l {z<£Xt} 
u({z : z < ext}) 



and finally Xt = v{{z : z < ext})- With a decomposition similar to (38)-(40), it is easy to show 
that the option price E[(X t — k t ) + \ admits an upper bound 

E[(X t - k t ) + ) = E[(X t - +Y t - k t )+] < E[(X* + j t t + Y t - k t ) + ] 

< E[(Y t + X t- fct) + ]P[X* < xt] + E[X t 2 ]3P[X« > Xt ]i 
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and a lower bound 

E[(X t - k t ) + ] = E[(Xt +Y t - k t )+] > e-~ Xtt n(Xl + itt + Y t -k t 
> e-~ Xtt ¥[Xt > -XtMi-Xt + Itt + Y t - k t )+] 



Similarly to (41)-(43), we have 

fk ■= 



(-extfi) 



z 2 v{dz) 



2 - a. 



<ext) 



2-a- 



and with the same logic as in ( 45 ) , we have that 



[*? > Xt] < 



1 



1\ 2e 



t ->• 0. 



\et, t tj 

It is now clear that one can choose e > so that the square root of this expression becomes 
equal to o(tk t a+ ). Since F[Xf < —Xt] admits the same estimate, and t\t — > as t — > 0, we get 
that 



E[(X t -h) + ]> m t E[(Y t - X t + Itt - h 



E[(X t - k t ) + ] < M t E[(Y t + X t- h) + ] + o{tk\- a+ ) 



where mt and Mt converge to 1 as t — > 0. Since xt = o(kt) and 7ft = o(kt), from (36), we then 
get 

+1 , ™ t{k t + Xt-ltt) l - a+ C + X _ a 



E[{X t - k t )+] > m t E[(aW t - X t + Itt - k t ) + ] + m t 



a+ — 1 



+ o(tk t 



m t E[(aW t - X t + Itt - k t ) + ] + -± + o{tk\~ a ' 

a+ — 1 

t(k t - X t) 1 - a+ c+ 
a+ — 1 

M t E[(aW t + X t- h) + ] + ^ ^ + o(tfc i-oH 

a + — 1 



E[{X t - k t ) + ] < M t E[(aW t + X t- h) + ] + M, 



Finally, since we also have x* = o(y/i) and 7ft = o(y / t), we get that E[(crPVi — xt + 7*i — ~ 
E[(aWt - k t ) + ] and E[(erWt + Xt ~ h) + ] ~ E[((rWi - &t) + ], which allows to complete the proof 
of Theorem [U 

5.2. Proof of Proposition [2} We first concentrate on the proof of Q. Let {a 2 ,u,b) be the 
characteristic triplet of X with respect to zero truncation function, meaning that 

X t = bt + aW t + J2 AX s> 

s<t 

where as usual for any s > we define AX S = X s — X s _ . 
Assume first that a = 0. The left-derivative of the function 



x i y (e 



-kt _ e x ) + 



IS 



x i y -e x l 



{x<-k t }i 
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and hence Ito-Tanaka formula [16^ Ch. IV, Thm. 70] applied to the process (e " — < 
(e 

J(o,t] 



kt _ y i e ids 



+ 



l(o,t] 



0<s<t 



e Xs -l 



ids + 



+ e Xs ~ 1 



{X s -<-k t }< 



0<s<t 



a X s -+AX s 



(e~ kt - e Xa 



for any t > 0, since, in this case, X has paths of finite variation. Since (AI s ) s >o is a Poisson 
point process with intensity measure v(dy) x ds, and X s ^ ^ X s for at most countably many 
time s in the interval [0, t] almost surely, taking expectations on both sides of the path-wise 
representation above and applying the compensation formula for point processes yields 

(47) E[(e-' 



-k t _ e x t)+] 



-ME 



e Xal {x s <~k t }ds 



+ E 



f'f i 

Jo Jr\{o} l 



< e -h _ e X s + V) + _ ( e -fct _ e A s) + j v ^ d& 



From Theorem 43.20 in [20], we have that 4* — > b almost surely as t — > 0. Therefore, for any 
£ > 0, each path X(u) satisfies the following inequalities 



Xt(co) > t(b — e) > — k t for all small enough t > 
(recall that by assumption k%jt — > oo as t \. 0). Furthermore, since kt J, as i J. 0, for all 



sufficiently small i we have X s (lo) > — for all s < t. Therefore it holds j f Q e Xa l{x s <-k t } 
almost surely. Since on the other hand we have 



ds -> 



e Xs l 



o 



{A s <-fe} 



e fct ds = e" 



the dominated convergence theorem implies 

ft 



E 



e Xs l 



o 



{X a <-k t } 



ds 



o(t) as i — > 0. 



To deal with the second term in (47), observe that for any e > 0, each path X(uj) satisfies the 
inequalities (b — e)t < XtiyS) < (b + e)t for all t sufficiently small. Therefore X(uS) also satisfies 
the following inequalities for any y £ R \ {0} and all sufficiently small times t > 0: 

( e -fct _ e ^s(w)+j/\ + _ (-g-fct _ e A s (w)-j+ < ^ e -fci _ e {b-e)t+y^+ _ (p-kt _ e (b-e)t\ + 

and 



(b+e)U 
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The second term in both sides of the above inequalities is in fact always zero for sufficiently 
small t. Therefore we get the following almost sure convergence: 



^K\{0} 



,X s +y\+ 



-kt _ e X s 



"I v(dy)ds 



(1 - e y ) + v{dy) as 1 1 0. 



is Lipschitz with a Lipschitz constant that does not 



Since the function y i— > (e~ kt — e Xs+y ) + 
depend on the path X(u), the dominated convergence theorem and the representation in (47) 
yield 

lim -E[{e- kt - e Xt ) + ) = [ (1 - e y ) + u{dy). 

40 t JR\{0} 

Assume now that a > 0. Define 

Z t := (b + y") t + AX s 
x " 0<s<t 



f(t,x) :=E 



1 - e 



x+kt+aWt-^-t 



and 



and note that 

(48) E[{e- kt - e Xt ) + ] = e~ kt E{f{t,Z t )}. 

The derivative of / with respect to x is given by 



(49) /'(*,*):= -E 
It can be computed explicitly as 

(50) f'(t,x) = -e x+kt N 



x+h+aWt-^t-, 

{x+kt+aWt-^tKO} 



> -1, 



5<T 2 t + X + fc 4 



(7 



v 7 * 



where N(-) denotes the standard normal CDF. Note also for future use that 



(51) f"(t,x) = -e x+kt N 



\a 2 t + x + k t 



1 

H i=n 



-\u 2 t + x + kt 



ay/I 



> -1, 



with n(x) = N'(x). 

Applying Ito's formula to the process f(t, Z) as a function of Z with t fixed, yields 

f(t,Z t ) = f(t,0)+ [ f'{t,Z 8 -)dZ.+ V [f(t,Z a )-f(t,Z a -)-f%Z.-)*Z, 
J(o,t] 



0<s<t 



f(t,o) 



+ ( b+ ^~) f t f'(t,Z a -)da+ V [/(t,V+AI s )-/(t,Z s _)] 

V 2 / ^ 0<s<t 



since AZ S = AX S for all s > 0. By taking the expectation and applying (48) we find 

ft 

-kt _Xt\+l _ _— fetipT-f /'-t n\l i „— fei 



(52) E[(e- fet - e At ) + ] = e" fct E[/(t, 0)] + e~ kt E 



f'(t, Z s )ds 



o 



+ e" fct E 



/' / {/(*, ^ + V) - f(t, Z s )} v{dy)ds 
Jo Jr\{o\ 
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The first term on the right-hand side of (52) is equal to the first term on the right-hand 
side of ([9]). As in the case a = 0, using the almost sure convergence — > b + the explicit 

00 as 1 \, 0, we get that 



form (50) of f'(t,x) and the assumption that — 



J / f'(t, Z s )ds -> 
1 Jo 



almost surely. Since \f'(t,Z s )\ < 1 for all t, s > by (49), the dominated convergence theorem 
yields 



E 



f(t, Z s )ds 



o(t). 



To treat the last term in (52), we use the fact that for any e > 0, each path Z(u>) satisfies the 
inequalities 

t(b + a 2 /2 - e) < Z t (u) < t(b + a 2 /2 + e). 
for all sufficiently small t. Therefore, since f"(t,x) > —1, the following inequalities hold 

(53) f'(t, t(b-e + a 2 /2) + dy) - 2te < f'(t, Z s (oj) + By) < f'(t, t(b + e + a 2 /2) + 6y) + 2te 

for any trajectory s h-> Z s (ui), where s £ [0,t], and all sufficiently small t. The random variable 
under the expectation in the last term on the right-hand side of ( 52 ) can be expressed as follows: 
r-t 



(54) 



{f{t,Z s + y)-f{t,Z s )}v{dy)ds = - 
iM\{0} 1 







ds / d9 







yf'(t,Z s + ey)u(dy). 



The path-wise bounds in (53) can be used to estimate (54) from above and below. For each 
path Z{uS) we have the following bound for y € (—00, 0) and all sufficiently small t: 



2te 



(-00,0) 



yv{dy)+ / d9 



! d6 [ yf'(t, t(b + e + a 2 12) + By)u(dy) 

JO J(-oo,0) 



1 



< 



< -2te 



ds / d6 







(-00,0) 



'(-00,0) 
yv{dy) + 



'(-00,0) 
yf'(t,Z t (u>) + 0v)u(dy) 



I dO [ yf(t, t(b-e + a 2 12) + 9y)u(dy). 

Jo J(-oo,0) 



The explicit form (50) of f'(t,x) implies that for all y < and 9 > we have 



f(t,t(b±e + a 2 /2)+0y) 



as t — > 0. 



Since f'(t,x) is bounded, the dominated convergence theorem yields 



d9 



(-00,0) 



yf(t,t(b±e+a 2 /2)+9y)u(dy) 



d9 



(-00,0) 



ye 9y v(dy) 



(-00,0) 



v(dy)(l- e y) 



as t J, 0. Formula (50 ) for f'(t, x) implies that for all y G (0, 00) and 9 > we have 

f(t,t(b±e + a 2 /2) + 9y) -> as t -> 0. 
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An analogous argument for y G (0,oo) to the one above and the representation in (54) imply 
the almost sure convergence 

1 « 



{f{t, Z s + y)- f(t, Z s )} v{dy)ds 



(1 - e y ) + u(dy) as t -> 0. 



t JO JR\{0} 

Finally, since f(t, x) is Lipschitz in x, with the Lipschitz constant independent of t, the dominated 
convergence theorem implies 



-E 



ft {f(t,Z s + y)-f(t,Z s )}v(dy)ds / (1 - e« 
Jo JR\{0} JR\{0} 



v{dy). 



This concludes the proof of ([9]). Note that in this proof, we did not use the condition in Q, 
but only the assumption J* R \r } \x\v(dx) < oo. 

We now concentrate on the proof of ^ . Since the Levy process X satisfies ^ , we can define 
the share measure IP, via §p\f t = e Xt , as in the proof of Theorem 3 Analogous to the equality 



in (30), we have 



(55) 



E[(e 



-k t _ g -AV 



where E denotes the expectation under the share measure P. Furthermore, it is well-known 
that under the measure P, the process X is again a Levy process with a characteristic triplet 
(<j 2 ,v, 7), where u(dx) = e x v(dx), and e~ x is a positive P-martingale started at one. The Levy 
measure v clearly satisfies 

|a;|z/(da;) < oo. 



Therefore we can apply ^ to the process —X under the measure P. Hence the identity in (55) 
yields: 

E[(e Xt - e kt ) + ] = e kt E[(e- kt - e aWt ~^) + ] + te kt [ (1 - e~ x )V{dx) + o(t) 

J(a,oo) 



E[( e ffWr *-V -e kt ) + ] +t 



(0,oo) 



l)v{dx) + o(t), 



where we used the Black-Scholes put-call symmetry given in d2§, the fact e kt = 1 + o(l) and 



the equality v(dx) = e x u(dx). This establishes the formula in Q and concludes the proof of 
Proposition [2} 

Appendix 

Lemma 5. Let X be a Levy process satisfying ^ and kt a deterministic function such that 



kt > Vt > and lim k t = as t J, 0. 

4 
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Then for any b G R we have 

E[(e Xt+bt - e fc >)+] = e kt E[{X t - k t ) + ] + 0(t) as 1 10, 
E[(e" fci - e Xt+bt ) + ] = e- kt E[(-k t - X t ) + ] + 0(t) as 1 10. 

Proof. Since < (Xt + bt — kt) + — (Xt — kt) + < 6 + i = 0(i), it is clearly sufficient to prove 
the formula for the call in the case 6 = 0. Let f(x, k) = (e x — e k ) + — e k (x — k) + and note the 
following: f' x (x, k) = (e x — e k ) + for all x G R and f' x \x, k) = e x ls x > k \ for all x G R \ {k}. By 
Taylor's formula we have f(x, k) = (x — k) 2 J^(l — 9)f"((l — 6)k + 9x)d9 for any x / k, and, 
considering fixed, we find 



E[f(X t ,kt)]=E 



(Xt-hf f\i-e)e k ^ Xt - kt h {kt+e{Xt _ kt) ^ kt} de\ <c Q nxh Xt \ 

Jo 



for some constant Co > 0. Under the assumption of the lemma, the right-hand side can be 
computed as 

mX 2 e Xt ] = -^E[e uXt ] 

OU z u=l 

A direct computation using the Levy-Khintchine formula then shows that ELY^e^*] = 0(t) as 
1 1 0. The put case is treated in a similar manner. □ 
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